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ABSTRACT 



A variational model for spatially dependent gamut mapping 
is described that includes inputting a gamut constraint, 
choosing a one dimensional gamut projection scheme, 
including selecting a transform color coordinate system, 
computing transform equations, and verifying gamut con- 
ditions. The model also includes inputting an original image 
to be rendered, where the original image is in a given color 
coordinate system, transforming the gamut constraint, the 
image, and the transform equations to the transform color 
coordinate system, whereby a three dimensional function is 
transformed into a one dimensional quadratic functional, 
finding a minimum solution to the functional, and trans- 
forming a projected image in the transform color coordinate 
system into a color coordinate system of a rendering device. 
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VARIATIONAL MODELS FOR SPATIALLY 
DEPENDENT GAMUT MAPPING 

TECHNICAL FIELD 
[0001] The technical field is color image recreation. 

BACKGROUND 

[0002] Gamut mapping is a method that maps a color 
image into a constrained color space of a given rendering 
medium, while attempting to preserve the visual quality of 
the image. The target color space is termed "the color gamut 
volume." Classical gamut mapping algorithms use fixed 
color mapping look up tables and ignore the spatial color 
configuration of the projected image. 

SUMMARY 

[0003] To improve color image rendering, a variational 
model for gamut mapping is described which does not 
require the target gamut to be convex. The model comprises 
choosing a one dimensional gamut projection scheme, 
including selecting a transform color coordinate system, 
computing transform equations, and verifying gamut con- 
ditions; based on a given gamut constraint; inputting an 
original image to be rendered, wherein the original image is 
in a given color coordinate system; transforming the gamut 
constraint, the image, and the transform equations to the 
transform color coordinate system, whereby a three dimen- 
sional functional is transformed into a one dimensional 
quadratic functional; finding a minimum solution to the 
functional; and transforming a projected image in the trans- 
form color coordinate system into a target color coordinate 
system. 

[0004] The halo elimination method may first be described 
in a one-dimensional formulation. Then, a true three-dimen- 
sional formulation can be reduced to the one -dimensional 
formulation. This reduction of the three-dimensional func- 
tional to a one-dimensional functional may also be used in 
other aspects of image processing besides halo elimination. 

[0005] In an embodiment, the one dimensional functional 
is applied to a continuous image representation. In an 
alternative embodiment, the one dimensional functional is 
applied to a discrete image representation. 

[0006] The variational model for gamut mapping may also 
be used to eliminate halos in a rendered image. Halos are 
artifacts that appear at object edges in the rendered image 
and are caused because there is no local restriction on the 
size of a projected image gradient. An embodiment of the 
variational model that solves the halo problem is based on 
the definition of a one dimensional gamut mapping problem 
in which the original image representation, denoted f(x), is 
assumed to be a positive, bounded and differentiable func- 
tion defined on a bounded interval. A projected image 
representation, denoted g(x), which is a positive and differ- 
entiable function, constrained by g(x)^c(x), is determined, 
with the objective that the projected image preserves the 
visual quality of the original image. This objective is trans- 
lated into the requirement that an average distance between 
the values of the two functions f(x) and g(x) as well as an 
average distance between their derivative values (represent- 
ing image edges) are small. To eliminate the halo phenom- 
ena, two additional constraints ((f) 2 ^g'f ^0) on the deriva- 



tive g' of g restrict the variations of the projected image to 
be smaller than the variations of the original image. This 
additional constraints agree with the intuitive notion of 
projecting an image into a smaller gamut volume. The 
additional constraints guarantee that at each pixel the deriva- 
tive of the projected image has the same sign but smaller 
value than the original image derivative. Given these con- 
straints, a one dimensional functional comprising an image 
distance term and a corresponding derivative-distance term 
can be solved to yield a projected image that eliminates the 
halo problem. Additional improvements, e.g., replacing f(x) 
with min(f(x),c(x)) in an image distance term and requiring 
that g(x)^min(c(x), f(x)) can be integrated into the one 
dimensional functional. Another embodiment of the varia- 
tional model that solves the halo problem uses a relative 
distance between the derivatives, rather than an absolute 
distance as in the prior embodiment. 

DESCRIPTION OF THE DRAWINGS 

[0007] The Detailed Description will refer to the following 
drawings, in which like numerals refer to like objects, and in 
which: 

[0008] FIG. 1 is a block diagram of an apparatus that uses 
a variational approach for gamut mapping; 

[0009] FIGS. 2A and 2B illustrate one dimensional image 
projections for non-convex gamut volumes; 

[0010] FIGS. 3A and 3B are a flowchart illustrating an 
operation of the apparatus of FIG. 1; and 

[0011] FIGS. 4A and 4B are a flowchart illustrating an 
alternate operation of the apparatus of FIG. 1. 

DETAILED DESCRIPTION 

[0012] Gamut mapping is used to modify a representation 
of a color image to fit into a constrained color space of a 
given rendering medium. For example, a laser-jet color 
printer that attempts to reproduce a color image on regular 
paper would have to map the photographed picture colors in 
a given color range, also known as the image "color gamut/' 
into the given printer/page color gamut. Conventional gamut 
mapping methods involve a pixel by pixel mapping (usually 
a pre-defined look-up table) and ignore the spatial color 
configuration. 

[0013] A variational approach for color image gamut 
mapping is based on an image difference measure which is 
a combination of two terms. The first term is the standard 
distance measure, aimed at preserving the original colors of 
the original image. The second term is the distance 
measure of the image gradient, which preserves the original 
image variations. The variational approach operates on each 
color channel, e.g., R, G and B, separately, and hence can be 
presented as a one dimensional (ID) approach. 

[0014] Considering one color channel only, an image S 
can be represented by a positive, bounded and differentiable 
function, f(xy), defined on a bounded domain Q. To project 
this representation into a restricted color gamut, which 
reduces to projecting into a line interval in the ID case, 
requires finding a positive and differentiable function g(xy) 
that minimizes the functional: 
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f (ft * (/ - g)) 2 dxdy + a f (V (A * (/ - *))) 2 rf*rfy, 



[0015] subject to C!^g(x,y)^C h , where Q and C h are 
constants that describe the 1 D target gamut restriction, 
k(xy) is a smoothing kernel, and a is a parameter. This 
minimization problem can be solved by an efficient qua- 
dratic programming (QP) iterative method. Also, and with- 
out loss of generality, the gamut constraint may be regarded 
as one sided such that g(x,y)^c. 

[0016] This variational approach to gamut mapping suf- 
fers from two main drawbacks. First, the efficient QP 
solution is based on separating the original three dimen- 
sional problem to three independent one dimensional prob- 
lems, with each of the one dimensional problems having a 
simple constraint of the form g(x,y)^c. However, real color 
gamut volumes differ significantly from a simple box in any 
of the common color spaces. Hence, a real gamut mapping 
problem cannot be solved as three independent one dimen- 
sional problems. In fact, gamut volumes are typically non- 
convex bodies and hence the a solution based on the 
variational approach may not converge. 

[0017] A second drawback with the variational approach is 
that the projected images produced by the variational 
approach have halo artifacts near image edges. Such halo 
artifacts are formed since there is no local restriction on the 
size of the projected image gradient. 

[0018] A corrected variational formulation for the gamut 
mapping problem eliminates the halo phenomena, and, 
given certain assumptions, reduces the three dimensional 
gamut mapping problem into a one dimensional QP prob- 
lem, even for non-convex color gamut volumes. 

[0019] For simplification, the discussion that follows con- 
siders scalar functions, defined on a bounded interval. An 
extension to color images, which are three dimensional 
functions defined on to dimensional domains, will be 
described later. 

[0020] An embodiment of the variational model that 
solves the halo problem is based on the following definition 
of the ID gamut mapping problem: Let the original image 
representation, denoted f(x), be a positive, bounded and 
differentiate function defined on a bounded interval. Find a 
projected image representation, denoted g(x), which is a 
positive and differentiable function, constrained by 
g(x)^c(x) with the objective that the projected image pre- 
serves the visual quality of the original image. This objective 
is translated into the requirement that an average distance 
between the values of the two functions as well as an 
average distance between their derivative values (represent- 
ing image edges) are small. In order to eliminate the halo 
phenomena, an additional constraint restricts the variations 
of the projected image to be smaller than the variations of 
the original image. This additional constraint agrees with the 
intuitive notion of projecting an image into a smaller gamut 
volume. The corrected variational formulation describing 
the gamut mapping problem now reads: Find the function 
g(x) that minimizes the functional: 

J(k*(f-g) 2 dx + aJ(Qc'<f-g))') 2 (2) 



[0021] subject to g(x)^c(x) and (l)21g'f ao. The second 
constraint ((f) 2 ^g'f i=0) guarantees that at each pixel the 
derivative of the projected image has the same sign but 
smaller value than the original image derivative. Additional 
improvements, e.g., replacing f(x) with min(f(x),c(x)) in the 
term and requiring that g(x) i min(c(x), f(x)) can be 
integrated into this formulation (i.e., equation 2). 

[0022] Another embodiment of the variational model that 
solves the halo problem uses a relative distance between the 
derivatives, rather than the absolute distance appearing in 
equation (2): Find the function g(x) that minimizes the 
functional: 



/(*. ( /-,))^ + a/(|g;-l)V 



[0023] subject to the constraints g(x) ^ min(c(x), f(x)) and 



[0024] This optimization problem can be simplified by 
defining a new variable 

""(**/)'' 

[0025] so that the above functional (equation (3)) 
becomes: 

J(k'(f-g)) 2 dx+aJ(u-l) 2 dx, (3a) 

[0026] and the set of constraints are: g(x) ^min(c(x), f(x)), 
(k*g)'-u(k*fHD and O^u^l. 

[0027] The optimization problems represented by equa- 
tions (2) and (3a) may be solved using standard discrete 
numeric schemes. The derivative operator in the continuous 
problems may be implemented in the discrete domain as a 
linear operation g'-Dg where D is a matrix. Hence, the 
discrete problems corresponding to equations (2) and (3a) 
are QP problems, i.e., the minimization objective functional 
can be represented as a quadratic expression in g and the 
constraints are linear in g. Thus, standard QP solvers can be 
applied for solving equations (2) and (3a). Explicit formu- 
lation of a gamut mapping problem as a discrete QP problem 
will be described later. 

[0028] Under certain conditions the three dimensional 
color gamut mapping problem corresponding to equation (2) 
can be mapped into a ID gamut mapping problem, even 
when the gamut volume is not convex. Considering the three 
dimensional formulation of the variational gamut mapping 
problem, the original color image may be represented in a 
three dimensional color space by three bounded and differ- 
entiable functions, denoted f 1 , r, f 3 , defined on a bounded 
two dimensional domain Q. Let V denote the restricted color 
gamut volume in this three dimensional color space. Three 
bounded and differentiable functions, denoted fp\ fj, 2 , fp 3 , 
are needed that minimize the functional: 
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[0029] subject to the constraints: 

(fp 1 . i/, fp^r. and (4a) 

WpW)so, (a/^ca/p'xa/^o for i=i, 

2, 3. (4b) 

[0030] In equations (4), (4a) and (4b), a' and, p 1 are 
parameters, k L are smoothing kernels, and 3^ d y denote 
partial derivatives with respect to the spatial coordinates x,y, 
respectively. 

[0031] The Euclidian metric used in equation (4) is con- 
sidered a good approximation to the human color distance 
measure only when special color coordinates such as Lab are 
used. 

[0032] The above-described functional (Equation (4)) may 
be used in imaging devices to map a gamut of an image to 
the gamut of a device that attempts to reproduce the image. 
FIG. 1 is a block diagram of an apparatus 10 thai may be 
used for gamut mapping according to the functional of 
Equation (4). An image capture device 11 receives an input 
original image and produces an output image S 12, which is 
an electrical signal representing a calorimetric value image. 
For example, the capture device 11 may convert at least 
three colorimetric values of each pixel of the original image 
into corresponding electrical signals. The electrical signals 
may indicate the L, a, b values, for example. Other colori- 
metric values may be the XYZ tristimilus value and the L, 
U, V or device dependent RGB values. A gamut mapper 15 
produces an in-gamut image S' 16. Finally, a rendering 
module 17 provides a rendered image S"18. The rendering 
module 17 may be implemented as a color laser printer. The 
thus-generated image S"18 may represent a best-fit image, 
given gamut limitations of the device in which the image 
reproduction is to occur. 

[0033] The image S 12 may represent the image as sensed 
by the capture module 13. The gamut mapper 15 applies an 
algorithm to extract and map the values of the image S 12 
into the gamut of the image reproduction device 17. In 
particular, the gamut mapper 15 may apply an algorithm that 
solves the problem represented by equation (4), thereby 
solving the gamut mapping problem and optimizing the 
output image S"18. 

[0034] The optimization problem defined by equation (4) 
is no longer a quadratic QP problem since, in the general 
case, the constraint represented by equation (4a) is not 
linear. Even if this constraint (equation (4a)) is approxi- 
mated by a set of linear inequalities, the gamut volume T 
may not be convex and then equation (4) may have local 
minima. To overcome the second difficulty, the projections 
into T may be restricted to one dimensional mappings. This 
restriction may produce a projected image that does not 
achieve the global optimum of equation (4). However, this 
restriction simplifies the solution to the gamut mapping 
problem significantly and performs better than classical 
gamut mapping algorithms, which are essentially one 
dimensional. 



[0035] The above-described ID projection scheme can be 
based on a gamut mapping projection method that satisfies 
the following conditions: 

[0036] 1. There exist three color coordinates denoted 1, 
h, s such that the original coordinates f 1 , f 2 , f 9 can be 
written in the form: 

fWCtOs+n'O^) for i-l, 2, 3, (5) 

[0037] where m £ and n { are differentiate functions. 

[0038] 2. The projection into T is performed along those 
lines in the this color space, where 1, h are kept constant 
and only the coordinate s is changed. The projected 
coordinates are thus lp=l, h p «h. An additional assump- 
tion of the variational model is that that s p ^S. 

[0039] 3. The gamut mapping volume T has the follow- 
ing property: For each point in the f 1 , f 2 , f 3 color space, 
that point's projection on the surface of T along the 
curves with 1, h constant, exists and is unique. There- 
fore, there exists a well defined continuous projection 
function s^s,!^) that maps the s coordinate of each 
color point outside of T to the s coordinate of a unique 
point on the surface of T (and for a point s, 1, h inside 
T: Sr<s,l,h)=s). 

[0040] One of ordinary skill in the art will readily appre- 
ciate that ID projections that satisfy the above conditions 
may be defined even for non convex gamut volumes. For 
example, if the original color space is the Lab space with s, 
h, 1 defined to be the standard spherical coordinates (where 
the L axis plays the roll of the z axis in the standard Cartesian 
system), then: L«s cos(l); a=s sin(l) cos(h); b=s sin(Z) 
sin(h) and the projection is along the radius to the origin, as 
illustrated in FIG. 2A. In another example, if the original 
color space is the Lab space and 1, s, h are defined to be the 
standard cylindrical coordinates (where the L axis plays the 
roll of the z axis in the standard Cartesian system), then: 
L=l; a=s cos(h); b=s sin(h) and the projection is performed 
along lines parallel to the (a,b) plane, as shown in FIG. 2B. 

[0041] As will now be described, the optimization prob- 
lem represented by equation (4) may be transformed into a 
ID QP problem when the projection from the image color 
coordinates f a , f 2 , f 3 to the new coordinates fp 1 ,^ 2 ,^ 3 
satisfies conditions (1), (2) and (3) defined above. In the new 
color coordinates hp=h, lp=l and hence the only unknown 
function is Sp(x,y). 

[0042] In the original color space the constraint repre- 
sented by equation 4(a) is (fp 1 , fp 2 » £p 3 )eT- ^ toe represen- 
tation of the original image in the new color coordinates is 
l(x,y), h(x,y), s(x,y), then the gamut volume constraint is: 
1 ami B p== 8 r(s,l,h). Thus, the unknown function s p (x,y) 
sliduld satisfy a constraint of the form: Sp(x,y)^c(x,y), 
where c(x,y)*s r <s(x,y), l(x,y), h(x,y)). This constraint is 
linear in s p (x,y). 

[0043] In the original color space the constraints repre- 
sented by equation 4(b) are: 

(^) 2 ^(a J / p i )(a J /)^o ( (a/)^(a/ p i )(V)^o for i i, 

[0044] Substituting equation (5) in the first constraint 
yields: 
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[0045] Denoting by fj.m^nj the functions 6^, d^dj) 1 , 
respectively, the previous inequality can be simplified to: 

[0046] which is a linear inequality in Sp. 

[0047] A similar computation, applied to the second 
inequality, which relates to the y derivative, yields: 

a y y-(r y i )Cn;)^(r y i )(m y vo' y f )( m t )a y 5 p ^-cr y t x«y i ). m 

[0048] The functional represented by equation (4) is qua- 
dratic in the original coordinates fj, 1 , fi, 2 , fp 3 and, according 
to equation (5), these coordinates can be written as a linear 
expression in s . Thus, the functional of equation (4) is 
quadratic in s p . The explicit expressions may be obtained by 
computing: 



(7a) 



(7b) 



f (k'tif-fpfdxdy* 
Jci 

\{k' * (m's - m's p )fdxdy = f(k' * (m'{s - s p ))fdxdy t and 

[{d x {k'*{f' -f' p ))) 2 dxdy = 

f \d x (k' *(m's-m's p ))) 2 dxdy = f (3 X (* ; *{m'{s - s p )))fdxdy. 



[0049] An operation 100 of the complete three dimen- 
sional gamut mapping algorithm to solve the optimization 
problem of equation (4) using the restricted ID projection 
will now be described with reference to FIGS. 3A and 3B. 
The operation 100 is a generic process for spatially depen- 
dent gamut mapping, and includes an off-line part (FIG. 
3(A) and an on-line part (FIG. 3(B)). Once the off-line part 
is complete, any color image may be processed according to 
the on-line part. 

[0050] The operation 100 begins in block 105. In block 
110, the gamut constraint is input. In block 115, a ID gamut 
projection scheme is determined by choosing the color 
coordinates 1, h, s, deriving m l (l,h), and n l (l,h) for i=l, 2, 3, 
from the transformation equations f^mX^bjs+n'fljh), and 
verifying that projecting along the s coordinate with respect 
to the given gamut volume T satisfies conditions (1), (2) and 
(3) noted above. Next, in block 120, the representation of the 
gamut volume is computed in the new coordinates and the 
ID projection function Sr(s,l,h) is derived. 

[0051] In block 125, the new representation of the deriva- 
tive constraints of equation (4b) is computed by substituting 
the transformation f^mXl.hJs+nfyh) for i=l, 2, 3 and using 
equations (6a) and (6b). In block 130, the 3D functional (4) 
is transformed to a ID quadratic functional in the function 
s p (x,y), by substituting the transformation f^m^ljhJs+nX^h) 
for i«l, 2, 3 and using equations (7a) and (7b). The new 
functional has the form: 



-continued 

pfaW SpWfdxdy} 



[0052] The operation 100 may now be considered com- 
plete in terms of establishing baseline conditions for pro- 
cessing a specific image S. The operation 100 then moves to 
block 135 for image processing of the specific image S. 

[0053] In block 135, the image S is input. In block 140, 
compute the representation of the original image in the new 
color coordinate l(x,y), h(x,y), s(x,y). In block 145, substi- 
tute the 1, h, s representation of the image S in the projection 
function Sr(s,l,h) and compute the constraint function c(x, 
y)-s r (s(x,y), l(x,y), h(x,y)). 

[0054] In block 150, find the function s p (x,y) that mini- 
mizes equation (8), subject to the gamut constraint Sp(x, 
y)=c(x,y) and the constraints of equations (6a) and (6b), 
using a standard QP solver. The optimal projected image is 
represented in the 1, h, s color space by lp(x,y)=l(x,y), 
h p (x,y)«h(x,y), s p (x,y)=s p (x,y). In block 155, transform the 
projected image from the representation lp, hp, s p back to the 
fp 1 ,^ 2 ,^ 3 color representation. In block 160, output the 
image S. In block 165, the operation 100 ends, 

[0055] FIGS. 4(A) and 4(B) illustrate an operation 200 in 
which the functional of equation (4) is solved using the three 
dimensional gamut algorithm based on the intensity pre- 
serving ID projection shown in FIG. 2B. The operation 200 
is a more specific example of the operation 100 shown in 
FIGS. 3(A) and 3(B), and relates to color transformation in 
the Lab color space. An original image S is represented in 
the Lab color space by (L(x,y), a(xy), b(x,y)) and the gamut 
volume surface, denoted 3r, is described in the same color 
space by a parametric representation (L r (u,v), a r (u,v), b r (u, 
v)). The operation 200 includes an off-line part (FIG. 4(a)) 
and an on-line part (FIG. 4(6)). Once the off-line part is 
complete, any color image in the Lab space may be received 
and processed according to the on-line part. 

[0056] In the operation 200, the image coordinates (x,y) 
are assumed to be discrete and hence the continuous for- 
mulation of the operation 100 is converted, when required, 
to a discrete representation. In this discrete representation, 
the derivative operators d x ,B y are replaced with standard 
matrix representations D x ,D y respectively and the integral 
formulation of equation (8) is substituted with a correspond- 
ing quadratic form. For a discrete image w the notation w x 
is an abbreviation for D x w and the operator of point-wise 
multiplication with w is denoted by (w). Using these nota- 
tions the gamut algorithm operates as follows: 

[0057] Referring to FIG. 4(A), the operation 200 starts in 
block 205. In block 210, the gamut constraint is input. In 
block 210, new color coordinates 1, h, s are chosen as: 



T(pf [d x {k'*{m'{s-s p )))?dxdy* 



(8) 



l = Uh~ arccosl 



(9a) 



[0058] and the inverse transformation is therefore: 

L-i; as ain(A); bs cos(/i). (9b) 
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[0059] The derived 1 D projection is performed along lines 
that are parallel to the (a,b) plane and intersect with the L 
axis. The gamut surface dT is assumed to have only one 
point of intersection, in each quadrant, with such lines. 

[0060] In block 220, the gamut restriction is transformed 
to the 1, h, s color space using equation (9a). The new 
representation is (l r (u,v), h r (u,v), s r (u,v)). The ID projec- 
tion function s r (s,l,h) is calculated by the following proce- 
dure: For each pair (l,h) find the coordinates (u(l,h), v(l,h)) 
that solve the system M r (u,v); h«h r (u,v). By condition (3), 
this system has a unique solution and hence there is a 
well-defined projection function Sr(s,l,h)=min(s, s r (u(I,h), 
v(l,h)». 

[0061] In block 225, the derivative constraints (6a) and 
(6b) are transformed, and in discrete notation read: 

( fly )^( fly )[(fc y )(oos(A)Msin(/i)A> P ^0 

^^^(^[-(^(^(A^Hcos^D^^O (6c) 

[0062] Finally, in block 230, direct computation of equa- 
tion (8), in discrete representation, yields: 

(s-s^CH.+HbXs-Sp) 

[0063] where H a =(sin(h))K B T (a a I+p a D JC T D x +|3 fl D y T . 
D y ' I )K a (sin(h)) 

[0064] and H^Cco^^a^+p^/D^P^/- 
D^cosCh)) 

[0065] 1^ and Kt, are the discrete representations of the 
convolution kernel k. 

[0066] FIG. 4(B) represents the on-line part of the opera- 
tion 200. In block 235, the original image S is received and 
in block 240, is transformed to the I, h, s color space using 
equation (9a). The coordinates of the new representation are 
denoted: l(x,y), h(x,y), s(x,y). 

[0067] In block 245, substituting the l(x,y), h(x,y), s(x,y) 
representation of the original image S in the projection 
function Sr(s,l,h) yields the constraint function c(x,y)- 
srCs^y), l(x,y), h(x,y)). 

[0068] In block 250, denoting sT(x,y) as the image that 
minimizes equation (8), subject to the gamut constraint 
Sp(x,y) ^ c(x,y) and the derivative constraints of equation 
(6c), the optimal projected image is represented in the 1, h, 
s color space by: lp(x,y)-l(x,y), hp(x,y)=h(x,y), s p (x,y)= 
s p (x,y). 

[0069] In block 255, the Lab representation of the pro- 
jected image is obtained using equation (9b): 

Lp-lp a p s p sin(/ip); £ p -J p cos(Ap) 

[0070] In block 260, the operation 200 ends. 



In the claims: 

1. A method for rendering color images using spatially 
dependent gamut mapping, comprising: 

inputting a gamut constraint; 

choosing a one dimensional gamut projection scheme, 
comprising: 



selecting a transform color coordinate system, 

computing transform equations, and 

verifying gamut conditions; 

inputting an original image to be rendered, wherein the 
original image is in a given color coordinate system; 

transforming the gamut constraint, the image, and the 
transform equations to the transform color coordinate 
system, whereby a three dimensional quadratic func- 
tional is transformed into a one dimensional quadratic 
functional; 

finding a minimum solution to the functional subject to 
the gamut constraint and a derivative constraint; and 

transforming a projected image in the transform color 
coordinate system into a color coordinate system of a 
rendering device. 

2. The method of claim 1, wherein the three dimensional 
functional is given by: 



fiftf**' *</' " f;)) 2 dxdy + f?£{d y k' •(/' - f' p )?dxdy} 

3. The method of claim 2, wherein the functional is 
subject to the constraints: 

(WW 1 . and 

(aJT*{a&) (a x f)so(a y £ i )g(3 y f p i xa x f)go for i-i, 2, 

3, and wherein a}, ft 1 are parameters, k are smoothing 
kernels, f is an original coordinate, f^ 1 is a projection of 
the original coordinate, and d x , d y denote partial deriva- 
tives with respect to the spatial coordinates x,y, respec- 
tively. 

4. The method of claim 1, wherein the transformed, one 
dimensional functional is given by: 



§ W/*' * imiS ~ Sp)))ldxdy ) + 

dxdy\ 



subject to: 

wherein Sp is a projection of s onto a target gamut and m* 
and n* are differentiable functions. 

5. The method of claim 1, wherein the color coordinates 
are 1, h, s, and wherein the transform equations are f=mXl, 
hJs+nXUh). 

6. A method for rendering color images using spatially 
dependent gamut mapping, comprising: 

receiving an original image representation, wherein the 
original image representation is defined by a positive, 
bounded and differentiable first function f; 
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mapping the original color image into a gamut space, 
wherein the mapped image representation is a positive 
differentiable second function g; and 

first minimizing a distance measure between the first 
function f and the second function g. 

7. The method of claim 6, further comprising second 
minimizing a distance measure between a derivative f , of 
the first function f and a derivative g'of the second function 
g, wherein (f 1 ) 2 ^g'f^O. 

8. The method of claim 7, wherein the first and the second 
minimizing steps are completed to minimize a one dimen- 
sional functional: 

JCk'(f-g)) 2 dx+a((k-(f-g))) 2 dx. 

9. The method of claim 6, further comprising second 
minimizing a relative distance measure between a derivative 
f of the first function f and a derivative g* of the second 
function, wherein g=min(c,f) and 



wherein k is a smoothing kernel. 

10. The method of claim 7, wherein the first and the 
second minimizing steps are completed to minimize a one 
dimensional functional: 



11. A method for rendering color images using spatially 
dependent gamut mapping, comprising: 

receiving an input color image in a Lab space; 

defining new color coordinates 1, h, s according to trans- 
form function 



12. The method of claim 11, wherein image coordinates 
(x,y) of the image are assumed to be discrete, and wherein 
the one dimensional functional is of the form 

where H^mCn))^ 1 " (aM+P'D/D^P^/D^K, 
(sin(h)) 

and Hb-CcosChJ^a^+p^D.+p^/D^Kfe (cos(h)) 

wherein D x , D y are matrix representations of derivatives, 
1Q, Kb are matrix representations of convolution ker- 
nels, and a\ p* are parameters. 

13. A computer readable medium, comprising program- 
ming for spatially dependent gamut mapping, wherein the 
mapping is carried out according to the steps of: 

receiving an original image representation, wherein the 
original image representation is defined by a positive, 
bounded and differentiable first function f; 

mapping the original color image into a gamut space, 
wherein the mapped image representation is a positive 
differentiable second function g; and 

first minimizing a distance measure between the first 
function f and the second function g. 

14. The computer readable medium of claim 13, wherein 
the mapping further comprises second minimizing a distance 
measure between a derivative f , of the first function f and a 
derivative g' of the second function g, wherein (f)2^g'f ^0. 

15. The computer readable medium of claim 14, wherein 
the first and the second minimizing steps are completed to 
minimize a one dimensional functional: 

JCk-tf-^^dx+ap^f-g)) 1 ) 2 ^. 

16. The computer readable medium of claim 13, wherein 
the mapping further comprises second minimizing a relative 
distance measure between a derivative f of the first function 
f and a derivative g' of the second function, wherein 
g=min(c,f) and 



inputting a gamut constraint for a gamut volume; 

transforming the gamut volume, the image, and the gamut 
constraint using the transform function, whereby a 
three dimensional functional is transformed into a one 
dimensional quadratic functional; 

finding a minimum solution to the one dimensional func- 
tional, wherein a projected image is derived; and 

transforming the projected image in the new color coor- 
dinate system into the Lab space color coordinate 
system. 



wherein k is a smoothing kernel. 

17. The computer readable medium of claim 16, wherein 
the first and the second minimizing steps are completed to 
minimize a one dimensional functional: 



/(*.(/- 8 ,,^ +0 /(^-i) 2 ^. 



* * * * * 
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